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THE  STUDY  OP  SUPEBSONIC  PLOW  ABOUND  DELTA 
WINGS  WITH  POBCED  ANTISYMHETBY  TAKING  INTO 
CONSIDEBATION  THE  PALLING  OFF  OF  FLOW  AT  THE 
LEADING  EDGES 


ST.  STAICU 

BUCHABEST  POLYTECHNIC  INSTITUTE 

This  paper  studies  the  supersonic  flow  around  thin  delta  > 
wings  haying  forced  antisymmetric  distribution  of 
incidences,  taking  into  consideration  the  separation  of 
flow  at  the  leading  edges.  Considering  an  Imaginary 
thin  wing,  equivalent  to  a real  one  from  an  aerodynamic 
point  of  view,  the  distribution  of  pressure  and  aero- 
dynamic characteristics  are  determined. 

1.  PBELIMINABY  CONSIDEBATIONS 

In  that  which  follows  we  will  do  a study  of  the  supersonic 
flow  around  thin  delta  wings  having  forced  antisymmetric 
distribution  of  incidences,  taking  into  consideration  the  falling 
off  of  flow  at  its  subsonic  leading  edges.  The  antisymmetric 
distribution  of  the  incidences  or  of  the  vertical  velocities 
corresponds  to  a formed  antisymmetric  curved  delta  wing  according 


to  * liniar  function  or  an  antisymmetric  deflection. 

While  the  incidences  are  very  small,  the  pressure  on  the 
wing  follows  approximately  the  stability  laws  in  the  hypothesis 
of  small  disturbances,  outside  the  system  in  immediate  proximity 
to  the  leading  edges,  where  the  corresponding  velocity  and 
pressure  are  finite,  though  infinite  values  would  result  from 
liniar  theory.  We  can  say  therefore,  that  outside  this  area  limited 
by  the  angle  of  the  leading  edge,  this  theory  is  valid  for  the 
whole  wing,  such  that  the  falling  off  of  flow  becomes  hardly 
felt  and  this  influence  is  small  on  the  whole  contents  of  the 
wing. 

As  the  values  of  the  size  of  the  local  Incidences  become 
greater,  the  flow  separates  at  the  leading  edge,  as  with  the  plane 
delta  wing  with  constant  incidence,  creating  a vortex  layer 
which  sits  above  and  below  the  wing,  producing  an  antisymmetric 
movement.  The  vortex  layer,  having  sufficiently  small  thickness, 
can  be  considered  a vortex  sheet  which  is  rolled  up  in  the  form 
of  a horn,  composed  of  a concentrated  nucleus  and  a marginal 
vortex  sheet,  starting  at  the  leading  edge. 

The  Incidence  being  variable  on  the  wing  surface,  the  axis 
on  which  the  horn  is  wrapped  will  be  a curve,  and  the  vortex 
generation  intensity  of  the  nucleus  is  variable  along  the  axis 
proportional  with  the  square  of  the  opening  of  the  wing.  For 


2. 


simplification.  In  that  which  follows,  the  axis  on  which  the 
nucleus  of  the  vortex  is  situated  is  considered  to  be  a 
straight  line.  Thus  the  system  of  two  concentrated  vortexes, 
of  the  same  intensity  and  sign,  situated  antisymmetrically  in 
reference  to  the  axis  of  symmetry  0X1  (fig.  1)  at  the  abscissa 
c and  the  ordinate  t,  will  bring  essential  modifications  on 
the  field  of  flow  around  the  wing. 


I'iB  1 


The  resulting  flow,  which  will  become  more  complicated, 
will  be  studied  on  the  basis  of  the  conical  theory  of  motion 


with  constant  and  antisymmetric  Incidence,  respectively,  {forced 
antisymmetry)  In  reference  to  the  axis  Ox*,  which  led  toward 

[ 

conical  motion  by  the  first  order  (conical  as  a matter  of  fact). 

We  will  allow  that  tht*  effect  of  the  falling  off  of  flow 
at  the  edge  of  the  wing  and  the  formation  of  the  two  antisymmetric 
nucleuses  consists  of  the  modification  of  the  field  of  vertical 
and  longitudinal  velocities,  having  as  a result  the  avoidance 
of  Infinite  velocities  at  the  leading  edge,  as  results  from  the 
classical  liniar  theory.  But  it  can  be  allowed  that  the  effect 
of  the  longitudinal  velocities  of  disturbance  can  be  substituted 
through  a corresponding  distribution  of  the  vertical  velocities. 

By  that  we  will  consider  a distribution  by  Incidences  or  by 
vertical  velocities,  forced  antisymmetrically,  so  as  to  correspond 
to  a real  case  of  an  Imaginary  thin  delta  wing  with  variable 
Incidence,  different  on  both  its  sides,  having  at  the  sane  time 
finite  velocities  at  the  leading  edge. 

It  will  be  allowed  that  the  real  thin  wing,  which  has  In 
a certain  way  finite  velocities  at  the  edge  through  the  effect 
of  the  separation  of  flow.  Is  equivalent  from  an  aerodynamic 
point  of  view  with  an  Imaginary  thin  wing,  having  the  same 
variation  of  Incidences  which  I defined  earlier. 

In  order  to  study  the  flow  more  easily  through  the  conical 

4. 


motion  method,  we  will  take  apart  the  imaginary  wing  corresponding 
to  distribution  of  the  above  vertical  velocities  into  three 
wing  components,  in  the  same  way  as  in  (2)  - (4)t 


1)  the  thin  wing,  having  a variation  of  forced  antisymmetric 
Incidence  suitably  chosen  In  order  to  follow  in  some  measure  the 
phenomena  of  the  modification  of  pressures  and  of  the  vertical 
velocities  of  disturbance  on  the  surface  of  the  wing  near  the 
leading  edge*  Thus  an  imaginary  thin  wing  is  obtained  with 
finite  velocity  at  the  leading  edge,  but  equal  and  of  opposed 
direction  on  the  two  sides,  higher  and  lower; 

2)  The  wing  of  •symmetrical"  thickness,  having  variable 
slope  in  the  same  way  as  the  Incidence  of  the  first  wing.  This 
wing,  combined  with  the  first,  will  form  a wing  with  different 
pressures  on  the  two  sides,  as  it  is  in  reality; 

3)  The  third  wing  will  have  symmetrical  "thickness",  with 
variable  slope  and  forced  antisy  :unetry,  however  in  such  a way 
that,  combined  with  the  wing  from  2),  a nought  mean  thickness 

Is  obtained,  characteristic  of  a real  thin  wing.  Superposltlonlng 
the  flow  around  the  three  wing  components,  we  will  obtain  the 
resulting  imaginary  wing,  equivalent  from  an  aerodynamic  point 
of  view  with  the  delta  wing  with  the  separation  of  the  flow  at 
the  edge. 


2.  THE  DETERMINATION  OF  THE  AXIS  OF  DISTURBANCE 
VELOCITIES 

5. 


In  continuation  we  will  follow  the  way  of  determining  the 


axis  of  disturbance  velocities  for  the  three  wing  components 
with  forced  antisymmetry,  being  necessary  for  the  determination 
of  the  distribution  of  pressures  and  of  the  aerodynamic 
characteristics  of  the  resulting  imaginary  wing,  which  are 
presupposed  to  be  the  same  as  real  thin  delta  wings,  having 
the  incidence  defined  by  the  relation 

w = 4;  ?r,„  r,  - 'f  :r,  al0J7«.  (1) 

Ue  will  note  further 

* . »:  = -«:  um,  = (2) 

the  vertical  velocities  and  the  incidences  respectively 

on  the  higher  surface,  w^,c?{  on  the  lower  of  the  thin  imaginary 
wing. 

The  movement  around  the  wing  being  conical  by  the  second 
order,  we  will  use  the  same  method  used,  considering  in  this 
direction  the  physical  plane  Oyz  (fig.  1)  normal  on  the  axis 
Oxj  and  having  the  coordinates 


the  axis  Oy  and  Oz  being  parallel,  with  0x2  and  Ox  3 respectively. 
Further  we  will  make  a similar  transformation  with  that  given 
by  Busemann  (fig.  2): 
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y . - - i*2(ua  + za)  , 

^ 1 _ $2  g%  & ~ ' 1 _ z't~  (X  — 'i  + >i)| 


l»i 


obtaining  a plane  which  has  the  property  of  keeping  the 
track  of  the  wing  (y  = y,  z *£■=*  0 ) in  the  true  magnitude. 
In  this  plane,  the  first  derivative  of  the  disturbance 
velocities  u,  v,  and  w are  harmonic  functions  and  can  be 
associated  to  corresponding  conjugated  functions,  in  such  a 
way  as  to  obtain  variable  complex  analytic  functions: 

* v f M- 

We  will  study  each  wing  defined  above  in  turn. 

2.1.  The  antisymmetric  thin  wing.  As  a result  of  the 
effects  of  the  two  nucleuses  of  the  vortexes,  the  vertical 
velocity  on  the  real  wing  is  modified,  as  well  as  the  first 
wing  component,  defined  above. 


Thus,  for  the  points  contained  between  ( -s<y<s  ) for 
the  track  of  the  wing  contained  in  the  plane  x « y + i^  (5)» 
the  vertical  velocity  will  be  considered  constant  for  xj  * constant 

7. 


w • 


• 


w'  = i J\  w'\l)  - (:  J-J  aio'  U<a,  (<;,* 

where  the  parameter  w|q)  corresponds  to  the  abscissa  y = s, 
and  for  the  area  y 6 [ -l,-sJuLe.\l  we  will  write 

w'  ~ ± — t j'i  a;0(i /)ir«,  (h 

\°ZJ 

such  that* the  leading  edge  (y  = ^1)  will  be  obtained 

M>U|  = ± jct  w,u  — f aio  P®  • *81 

This  continuing  variation  of  vertical  velocities  ( or, 
more  precisely,  of  the  parameter  w*g(y))  corresponds  to  the 
continuing  distributions  of  elementary  edges,  situated  on 
the  wing  In  the  Interior  of  the  considered  space,  which  gives 
each  point  y * TJ  the  elementary  drop. 

However,  taking  into  consideration  previous  works  M. 
&}.  the  contribution  of  elementary  edges  situated  In  the 
point  y * In  the  expression  of  the  axis  of  disturbance 
velocity  In  the  point  x,  applying  the  similar  hydrodynamic 
method,  it  will  be 


"“'if’  r;': 


These  contributions  of  the  edges  in  the  expressions  of 
the  axis  of  disturbaneeaaitd  vertical  velocities  on  the  wing 
are  realised  placing  on  the  trace  of  the  wing  from  the  plane 


x * y ♦ 1>  (5)  some  singularities  (sources)  by  order  two. 
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In  this  way  we  can  get  In  calculation  the  effect  of 
the  two  nucleuses  of  the  vortexes  concentrated  among  the 
continuing  distribution  of  sources.  This  division  will  be 
necessary  to  satisfy  the  conditions  Imposed  by  the  problem 
of  obtaining  concomitant  axis  of  disturbance  and  vertical 
velocities  Indicated  above  on  the  basis  of  observations 
and  the  resulting  experiments.  From  that  we  will  chose 
a llnlar  variation  of  Intensities  of  sources* 


9'(.V)  = q 


'{'-I') 


(»<>)<  l), 


which.  In  the  case  of  homogenous  motion  of  the  second  order 
(n  ■ 2) , Is  written 


7m  - ± ?’n|l  — ^ j>  q’2l  = ± j- 


The  axis  of  disturbance  velocity  for  the  thin  wing 
component  will  be  obtained  through  the  addition  of  the 
contributions  of  all  the  elementary  distribution  edges  and 
of  the  subsonic  edges  Cil  under  the  form 


An  f — Ql,  r (>os/( 

/ 2 i»2  rr 


+ - / )<«*  + dT1  " 


- l\[l  ~ 1 )(?’n  - Tq*n) ilr>' 


which,  after  the  accomDliahment  of  the  ealeulatl ona.  becomes 


I cos h 


*x 


U*  -f  r) 


— (?s* o - q\\*)  (»  + -r)^l  — 

+ — 29*,/)  cos" 1 * 1 9*i^|/  1 - *'  1 “ j’  J - ( I-1) 


In  whlch/U^Qi,, ^ , yxl  are  some  constants  which  will  be 

* 

determined  below, beln8  the  Intensity  of  a source  sitting 
In  the  origin  (x  * 0),  due  to  the  central  edges  which  appear 
at  the  wing  with  forced  antisymmetry. 

2.2.  The  wing  of  symmetrical  thickness  having  equal 
slope  with  the  incidence  of  the  first  thin  wing.  Through 
the  Introduction  of  this  wing  of  symmetrical  thickness,  the 
accentuated  peaks  of  pressure  on  the  lower  side  of  the  wing 
1 8 removed,  where  the  distribution  of  pressure  obtained  through 
the  superpositioning  with  the  first  wtng  component  will  be 
different  from  that  on  the  higher  side.  Following  the  general 
method  of  conical  motion  Cil.  for  a wing  of  symmetrical 
thickness  with  the  variation  of  slope  given  by  the  same 
distribution  of  sources  (11)  we  will  write  the  following 
expression  for  the  axis  of  disturbance  velocity: 


V/,  - fUy 


l if'  I")**  1 X«‘‘)COS'1  ‘fM 
i !(■  - ? H xr~dr‘ 1 


h ^ </:*  J-  cos  /< ' 1 y ± L. 


(Ill 


Accomplishing  the  Integrals  above,  we  will  find 

1 v/, 

^ 1 li(I  - a) 


l {<?*»  t '/*,  ' ) - l + ^006 


os  It  1 


1 
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I* 


- (*  - J-)  I 


l/l 

I 1 il  J B(»  ~ X)  J 

- (ti»  - </>)[(!  | x)(l  - L~  f]C(K  ,ri  J_+  _ 

L \ -l  J B(l  + x) 

~ (*  + •«  ) f l - s “ ?)  cos  A-I 

l a/  i /*(«  + X)  J 

— — 2{i 1 (s,t>  */y/  — sit 


s*«  'Bn)  (ql  - 2£,l)  - 


V M-** - n— SV)]}  + f eir^a-f ±x,  |15, 

where  L represents  the  contribution  of  the  subsonic  edge 
haying  the  slope  equal  wlthtt^X,  , and  the  term  which  presents 
the  coefficient—  appears  only  ln-the  case  of  forced 
antisymmetry. 

2.3.  The  wing  of  symmetrical  thickness  compensates 
for  slope.  The  Introduction  of  the  effect  of  the  wing  from 
point  2.2  makes  the  resulting  wing  have  "symmetrical  thickness". 
In  order  to  compensate  this  work,  we  will  Introduce  a new 
distribution  source  of  a certain  form,  which  will  restore  the 
wing  to  a mean  nought  thickness.  The  variation  of  the  vertical 
velocities  w*  given  by  these  sources  will  correspond  to  a 
"compensating  wing  slope",  having  still  a symmetrical  thickness 
In  reference  with  the  axis  0x3,  and  antisymmetric  face  for 
0x2.  This  wing,  having  at  the  edge  of  the  wing  the  velocity 
“i » w111  cancel  the  mean  Blppe  and  the  effect  L of  the 
slope  of  the  edge  of  the  wing  2.2.  The  distribution  of  sources 
of  Intensity  q*  will  be  necessary  to  create  on  the  lower  side 

11. 


J 


of  the  wing  a distribution  of  pressure  without  accentuated 
peaks,  approximately  llnlar,  with  the  exception  of  the  points 
near  the  leading  edges. 


To  simplify,  we  will  take  the  following  expression  of 
the  distribution  of  the  intensities  of  the  sources: 

(- l <'»<<)■  ('«) 

We  thus  obtained  two  “large  wings"  In  order  to  form 
a single  one,  havfcng  the  slope  variable  In  such  a way  for 
the  mean  to  be  nought. 


The  axis  of  disturbance  velocity  Ulc  for  the  third  wing 
will  be  the  following: 


(17) 


In  the  course  of  the  calculations  we  are  led  to  the 
expression  ~ at.  = ~ + ««  ;r)  cosh~ ‘ ^ ~ 


— (qV«  - <r)  <’Osli 


,H-  7?2  lx 

B[l  + *) 


+ 


+ 2B/|(?;;»in- '}«  l-  (1  - V t - 1 - K*.r*  + 

»■ -'l^l  + ~ g™'””  1 f 1 


(1«) 


nr 


^ »*» 

where  the  term £Tyv|  corresponds  only  to  forced  antisymmetry. 

Superposltlonlng  the  three  wing  components,  the  resulting 
Imaginary  wing  Is  obtained,  equivalent  from  an  aerodynamic 
point  of  view  with  the  real  wing,  for  which  the  axis  of 
disturbance  velocity  has  the  expression 

(Ux  = (Uu  + fUu  + (19) 

which  will  be  the  anti  symmetrical  face  of  the  axis  of  symmetry 
Oxj  and  continues  In  the  origin  0.  We  will  observe  that  the 
velocity  Uji  on  the  higher  surface  Is  equal  and  of  the  opposed 
sign  with  that  of  the  lower  surface,  as  corresponds  to  the 
thin  wing. 


3.  DETERMINATION  OF  THE  CONSTANTS 

For  the  determination  of  the  constants  qgo  and  q^i  which 
appear  In  the  expression  (13)  of  U11#  we  will  start  from  the 
following  conditions t 


lmC$ 

- ( 
(lx2  \ 

drUxl\ 

. dr)  J 

II 

(20) 

ImC(  1' 

d2 

/ d 

\ dx  = 0, 

(20') 

dx2 

{ dr) 

J 

deduced  from  the  theory  of  conical  motion  [l},  the  integration 
being  made  on  • semicircle  of  very  small  radius  around  a certain 


13 


point  y on  the  wing,  contained  in  the  interval  (y«a,  y-1) 
(fig.  2)  • 


Thus  we  will  obtain  the  relations 

, _ d<  r)  (2  - BV)  ' = I 

9ia  ~ dr]  (1  - Btr,iYl‘  ' dr)  (1  - BV)’'* 


(2i> 


which  stabilise  the  dependence  from  among  the  intensity  of  the 
sources  and  the  variation  of  the  vertical  velocities  on  the 
thin  imaginary  wing. 

Starting  from  these  relations  and  keeping  in  mind  (111, 
we  will  put  the  conditions  at  the  limit  in  the  points  7|  * s 
and  T|  ■ 1 for  the  parameter  Wj0  of  the  vertical  velocityx 

„ii»  _ „,!«»  = jL  j'  (<->?)(i-  W 


tPii,1  - Wtf  = -_s-  \ — . ; - - dr), 

I ).  rj  (2  - B V) 


(22) 


toil1  - C (l  -ft)  (l-  BVf*  dr), 

I 


(22') 


from  where  we  find  the  first  relations  among  the  constants 

tto,  • 


B*fi-  yi -B***-  COR*'1—  + COR  h~l  — 
2 Jil  Us 

— — ~ j Vl  — BV  )-  -jjj  (sin-1  Bl  — Rin-1  B»)  — 


/V. 


< 


r ' 


JpMi + . '+ m _ - « pM»  _ 

2 \ |/2  — m fi  4 fl(  - b* 


Kin 


-1 


idMU-1-  U-'  *-=-? 

^2  4-  B*  / ^2  Bl  1 \'i  - 


B)|(2  _ 


Bl 


. 1 + Bl\2  • 1 - M2  . . 14  \1 

-Bln  1 -rr= — Rin  1 — r 4 Bin  1 l7-=  — 


\2  4 Bt  !‘Z  - li * 

= trio  - «C, 


V‘2  i B» 


(23) 


7« 

8B 


^ f.3  ^sin-'TH  - sin"1®*  4 Bl  Vi  - B*i*  - /?#  Vl  - Bvj  4 
2 ( B/  (1  — B*(*)*'>-  Bs(1  - /?***)•'•)  4 -|-f(l  - B4(*)*'» 


- (1  - BV)*'» 


)] 


— «i!?> 


Wio  - ^io 


ID 


(23') 


However,  for  the  other  parts,  the  wean  vertical  vleoclty 

w xi  or  the  mean  Incidence  ^aXj  of  the  real  wing,  equal  with 
10 

that  of  the  first  wing  components  or  of  the  resulting 
Imaginary  wings,  is  obtained  takbng  Into  consideration  the 
two  large  wings  2)  and  3)  compensating  reciprocally  creating 
a mean  nought  elope  and  we  will  be  able  to  write  the  relation 

4-U.OD^vw  (21) 

I Jo 

only  for  the  thin  wing  1.1  Accomplishing  the  Integral, 
we  obtain 


M’>«  - toil,’  - (I  [ li  (til  — 5.)  (l  - BVf-  - 


15. 


-T 


1 


-IU(\-  - a|«/(i  - - *W*|- 

— m\  1 — h*l*  - /<«V i - v***  -Mu  ‘M~  siu  ‘i<«j|* 

ft 

For  the  determination  of  the  constants  A21  and  Q21 

which  appear  in  expression  (13)  of  the  axis  Telocity,  the 

✓ 

variation  of  Tertical  Telocity  w^qX^  will  be  taken  into 
consideration  from  a point  on  the  wing  to  the  one  of  nought 
vertical  Telocity  ( for  exaaple  on  the  Mach  circle).  As  in 
previous  work  £ 2 J - 07,  in  order  to  avoid  the  difficult 
calculations  which  appear,  we  will  consider  that  the  llnlar 
distributed  sources  in  the  interval  (s,l)  are  concentrated  in 
y > s,  for  intensities  Q20  and  Q211  in  such  a manner  that  we 


have 


a **  j’  - y-J- 


Proceeding  in  this  way,  we  will  write  the  relations 


•cercul  Midi  1 / I _ lit,.*  .12  //' 

-fail  U l J o, 

•'artpA  X lix* 


f t - 

iIj* 


where  is  the  axis  of  disturbance  velocity  for  the  simplified 

case  of  sources  concentrated  in  y * %%  given  by  the  expression 

««  - — <U\  - f !il-  + - (Qta  + V*.  -0  cos  h~l  {-—]  - 

j-,  H2  — x*  7t  r ai  («'  — x) 


— (V20  - Qm*)  oos/i  1 

TC 


Through  the  accompli shaent  of  the  integrals  which  appear 
above  on  the  axis  of the  ordinates  (y  ■ x a l^)  between  the 
Units  0 andoo  for  (2?)  and  on  the  axis  of  the  abscissas 
(y  - 0)  between  the  limits  1 and  l/B  for  (2?f,  result 

„ (2  — ti*  l*)  E (*)  — /?WT  (*)  9 . ( 

* I*  (1  -«*!*)  + 4«  f / — 

- li*s’*)E(k)  - HV*(K(k)  - Il(g,*))]  I V*.(j  - f)'  X 

X T^r£[(2  “ /?V*)  11  (p’  *>  - 7r(,,)l  + -Qhy  X 

1 — /f2S  2 L J 7T  I 


X 


-2/?(*)]=o, 


(29) 


O'  j _ 

2^9  -/{*((  f2*)2 


[ /?2/2  f 

i 2 -7-I  f 

1-  M%*n  4 3/f 

L l 

Ml 

w l 

(29') 


in  which  K(k),  E(k),  IT(p,k)  represent  the  complete  eliptioal 
integrals  for  the  first,  second  and  third  Instances  respectively, 
having  the  module  k and  the  parameter  p given  bp  the  relations 

ni(,D*HDt  [f  - *<»>*<*-*•>  + <*«*>  - 

-]?(*))  F(9o,  *')]•  (:i°) 

k = yr-^Ti'V*,  P = 7?V*  - 1,  *'  - /?/,  9o  sin  1 *-■  (30  ) 

17. 
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Due  to  the  separation  of  flow  at  the  edge  of  the  wing 
and  as  a result  of  the  presence  of  rortexes  on  the  higher  slde( 
finite  velocities  are  achieved  In  those  points.  Imposing 
this  condition,  we  will  be  able  to  write  the  relation 

A21  - 0.  (3D 

Eliminating  mC,  «>'»«  y;,  among  the  equations  ( 23)  • (23^. 

(25),  ( 29)  • (2/)  and  talcing  Into  consideration  (31).  we 
obtain  the  constants  <120  q-ji  : 

/...  • i*« W’io  1 (32) 

,;>  - - -f  -r:,,!  ,n  _ /„J..  - 

i 20  **  21  J 21  ° 20 

where  we  made  the  following  notations* 


l,„  = 111  1 2 |yi  - «»/•  - ^ - j j V 1 - - cos  * 1 ^ + 

.1.  tos/r1— 1 -1-  -i- [*2  (sin-1  HI  — sill-1  Bn)  — (HI  — ]/'i)  X 
a x | in  L 


( .1  -Blfl  1-B*1(2\ 

x sin-'  — , -- *»*  -./=  ,TT 


V2  - HI 


\'l  - Bn 


(K!  + V2)(»i.r'  <* 
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•'*  ^iu"1  Bl  — siir*  ft*  | Bl\  1 — ft*/*  - Bs\  1 — ft***  j + 


4 


\l[Bl(l  - if*/*)*'.  - B*(\  — ft2**)1'.]  -f 


5 ft/ 


'Mj 


8 (ft*/2/T(*)- 2E{k))  - ft2** 


- ft**'*)  II  ( p,  k)  - K(k)]  4 2E(k)  - ft*/*  Jf(*) j 

+ — - r»»)u  - fl’*1)5'2  - «/<l  - ft*/2)5'1  1 

- 2 [ ft/(  1 - ft2/*)3'2  - ft*(1  - ft***)*'*  ] 


- $(Mfi  — ft*/*  /?* fi  : ft*** 


sm 


ft/— sin1  ft*)]  • [XV") 


**  ** 

The  constants  q22  and  ^21  ®PPe*r  In  the  expression 

Ulc*  given  by  (18),  are  determined  taking  into  considemation 
the  rolft  of  the  third  wing  component,  which  will  hare  the 
mean  slope*  - wiO*i«  Similarly  with  (24),  we  will  write 

1 f* 

--I  tj)  «1  »j  — (•*•> 

* •'O 

19. 
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(3b') 


In  which  wio^  *8  the  ▼wtical  velocity  at  the  middle  of 
the  wing  compensating  for  slppe  which  is  obtained  starting 
from  the  relations  (16)  and  ( 21 ) * 


- i- 


= [1(1  - n2i2fn  I -'-(l 

1 f sin"'  It  l \ 2 I 

i (v*  - ”’*■ +~i<rr  r.r 


(JO) 


+ 


V\ 


4.  THE  DISTRIBUTION  OP  PRESSURE  AND  AERODYNAMIC 
CHARACTERISTICS 
21* 


i j— — 


We  have  shown  above  that  the  axis  of  disturbance 
velocity  of  the  real  wing  results  through  the  superpositioning 
of  three  imaginary  wing  components,  obtaining  formula  (19)* 

For  the  calculation  of  the  distribution  of  pressure,  the  total 
axis  velocity  given  by  (19)  (fig*  3)  will  be  considered: 


c = _ 2 — 1 
V 


Moving  along  to  the  calculation  of  the  coefficient  of 
lift  of  the  wing,  we  will  make  the  observation  that  the  wings 
of  symmetrical  thickness  do  not  give  liflt,  so  that  only  the 
coefficient  of  lift  given  by  "the  wing  lift"  will  be  taken 
into  consideration: 





Taking  into  consideration  (13)»  we  will  obtain  the  following 
expressions  ofthe  coefficient  of  lift  ofj^the  wing: 


+*l(’ -1- ?“"t- 

-tF;S 


S2  I _ . l 
4-  — cos  h 1 — 
l2  * 


3- 


The  coefficient  of  the  moment  is 


‘w“,  hi' “ £ {I vu + *[(!( l-?)- 


] au-io° 


Moo  =1.9 
Bl  -0.67 0 


»f  fifetsure. 

dm  mum  u msiuNi 


08  y/L 


where  H Is  a length  of  reference,  which  can  be  taken  equal 
with  the  unit. 


As  concerns  the  parameter  - which  enters  In  the  above 

1 

expressions  and  which  determines  the  limits  of  distribution 
of  sources,  we  will  observe  that  It  depends  on  the  position 


* > j 


of  the  center  2 of  the  nucleus  of  the  vortex.  In  previous 
works  (2),  (3)  we  analized  in  detail  the  problem  of  the 
position  of  the  nucleus  of  the  vortex.  In  this  manner,  taking 
Into  consideration  the  experimental  results  obtained  by 
various  authors  on  the  plane  delta  wing  with  constant 
Incidence,  we  will  allow  for  the  position  of  the  nucleus  of 
the  vortex  the  following  approximate  formula  of  variation 
with  the  Incidence: 

c 1 (44) 

i ~ i + i,i  <*r 

in  which  <x  represents  the  Incidence  of  the  wing  in  the  point 
which  would  represent  the  center  of  gravity  of  the  aerodynamic 
df fects: 

x,  » l . (45) 

L 4 

As  has  proceeded  and  in  previous  works,  we  will  consider 
the  following  relation  between  c and  s: 

c = * + 1(1-*);  =1,2^0.).  (46) 

5.  THE  PLANE  DELTA  WINO  WITH  FORCED  ANTISYMMETRY 

In  this  case  we  will  anallze  briefly  the  flow  around 
thin  delta  wings  with  the  Incidence  oc  equal  and  of  opposed 
sign  on  Its  two  halves,  taking  Into  consideration  the 

24. 
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the  separation  of  flow  at  the  edges. 

The  formation  of  vortexes  has  the  direct  result  of 
producing  a complex  field  of  vertical  velocities,  which 
will  modify  the  flow  in  such  a way  that  the  pressures  will 
be  finite  at  the  edges. 


Proceeding  as  above,  we  will  chose  for  the  distribution 
of  sources  the  following  llnlar  functions: 

Po-J  (-S<T)<s), 


<7'<ri) 


(•7) 


rje  [~  h — si  u [s,  l], 


In  the  case  of  the  first  two  Imaginary  wing  components,  and 

*"(*))  = *.]■  < — * < ti  e *),  M*) 

for  the  wing  of  symmetrical  thickness  compensating  for  slope. 


ij 


Applying  the  stability  formulas  In  the  exercise  of  conical 
motion  W.  the  axis  of  disturbance  velocities  of  the  three 
Imaginary  wings  will  be  the  following: 


-«*[(. -*)[l  _*-+-!) 
- **  7T  L \ 2 1 ) 


a x 
■ *) 


- («  + •») 


i2  + 


cos  h 1 — + a: 

l(v  + x) 


F:  ? 


COS 
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• v 


, vo  r *a  - **  ( cos  *-»  - cos h->  - + 

+ -[—  2T“l  !(•-*>  «(*+*)i 


+ * 


|/l  - ■ 


(19) 


for  the  thin  wing,  and 


fl„r  t i + *\  i — bhx 

‘ tt  LC  1 2<  / B(l-x) 

( l — x\  1 + B2lx 

— ll  -1-  x)  1 1 I cos  h 1 

1 h H 2 1 ) B(l  + x) 

( s + X\  1 — , 

_ (S  _ j-)  1 1 cos  h 1 h 

1 ’{  ‘21  ) B(s-x) 

8 -X\  1 + B2SX 

+ l ' '1  2i  ) B(8  + x) 

— i-  (sin-1  Bl  - sin-1  Bs)  xfi ~~ 1 + 

Bl  J 

n f v2  — X2  / 1 — B28X\  , 1 + B2SX\ 

« [ 2i  v B(8  - X)  I B(8  + X)  ) 

+ JLl IT~inj*x2  sin'1  7Js]  ±L, 

Bl  J 


(50) 


in  the  case  of  the  wing  of  symmetrical  thickness.  For  the 
third  wing  component,  "compensating  for  slope",  we  will  have 
the  velocity 


(U  = k°  - f B2  (l2  - x2)  ( cos  A 1 

* 2n  B2l  L V 


1 _ B*lx 
B(l  — x) 


cos  A-1 


1 + B2lx\ 

B (l  + x)  J 


+ 


-t-  2B  jfl  — B2x2  sin"1 


(51) 
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we  will  determine  the  constants  qQ,  p0,  k(j  in  a similar 
way  as  has  proceeded  and  in  previous  work  £2}  - C*7.  Taking 
into  consideration  the  relation 

o'  ^ W-  Aw'  . (r.2) 

«1y] 

known  from  the  theory  of  conical  motiony  and  by  the  relations 
(47),  (48),  we  will  put  the  conditions  at  the  limit  for  the 
the  vertical  velocity  of  the  first  wing  component: 

w0  - w = P—  C Kl  - B*  7)*  d7) , (53) 


where  wQ  and  Wj^  represent  the  vertical  velocities  at  the 
middle,  the  respective  edge  of  the  wing. 

Defining  a mean  vertical  velocity  equal  with  that 
corresponding  to  the  incidence  of  the  real  wing,  we  will  write 
for  the  thin  wing  the  relation 

tr  = --  (to' 7]  |"  + ir'rjj  — ( r)dw'  Tjdtr'V  (5-1) 

The  constant  A^,  which  we  determine  taking  into  consideration 
the  variation  of  the  vertical  velocity  from  a point  on  the 
wing  to  the  one  of  the  nought  vertical  velocity  (the  Hach  circle). 


2? 


as  in  previous  work  £2] , 03  * we  will  cancel  It,  imposing 
in  this  way  the  condition  of  nought  velocity  at  the  edges. 
By  accomplishing  the  Integral  remains 

i.cercul  Mach  tl'H.oo)  1/  1 L12  v» 2 

Im(  ^ — d-Kl-tr,,  (»5) 

Jarlni  (0.  I)  * 


where  is  the  axis  of  disturbance  velocity  in  the  simplified 
case  of  a concentrated  source  in  y ■ s'on  the  thin  wing. 

Through  the  accomplishment  of  the  integral  (55)  on  the 
axis  of  the  abscissa  between  the  limits  1 and  1/B  results 

~ -£)].  m 

in  which  s' and  Qq  are  given  by  the  relations 

' 3 9.(l-T)’  + P"5 

v„«“l (57,) 

Similar  with  (5*0*  we  will  have  in  the  case  of  the  wing 
compensating  for  slope 


---H”’"’1  [-[■“"■’I’ 


Integrating  the  equations  (53)t  (53l.  (5*0  and  (58), 
we  will  obtain  an  algebraic  system  of  equations  from  which 
the  constants  q0,p0,  kQ  will  result: 
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[(i-.Ms + &**]* 

"-  = l(yTzrhv*-]fi  _*'*  ) — LL — if — * H 

'•  r 1 " j(- 


— - if2/2  - hs  \ L - if2#2  + sin-1  Iil  - sin"1  if#  + 

+ .—((I -Bn*)** -a- b* »*)**)+  |-^(i  - (i  - if***)*'4)! 
3#/  3#/  gu  j 


+ 3 «t(" 


+ sin-1  if/— sin-1  if*  | — ((1  - if2/2)52- 


- (l  - if*#2)5'*)  + hi  J cos  irl  .y.  - cos  /«-•  ~JJ  if#  y 1 


+ sin-1  if# --(1  - (1  - if2#2)5'2)  l"1. 

3 hi 


if***  + 


Jt„  ^ _ gj(l  _ m*fV  - (1  - if2#2)5'2  -t-  ~ hi  \ hi  (l  - /f2f2),/2- 

l 2 


_ /}S(1  _ fl*#2)i*  + sin*1  iff  — sin-1  Bs]  + 


4-  - [1  - (i  - 

In 


tj  Cl  - (1  - 


if2/2)5'2! 


The  coefficient  of  lift  is  given  only  by  the  first 
wing  component: 


1 


-o-snMMijKr-i 

and  the  coefficient  of  the  moment  will  be 

' ~sk  i "• » - 


(01) 


Concerning  the  position  of  the  nucleus  of  the  vortex 
In  the  case  of  plane  delta  wings  with  forced  antisymmetry, 
the  Same  formula  Is  taken  as  In  previous  works* 


/ 1+1,7  («)»■■  ’ 


1,2 ----0,2. 
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